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ABSTRACT 

The  mechanical  properties  of  a  thin  layer  at  the  surface  of  aluminum 
alloys  change  progressively  during  fatigue.  The  alterations  are  more 
pronounced  and  rapid  for  fatigue  In  moist  than  In  dry  air.  We  characterize 
the  alterations  observed  as  a  change  In  the  "ductility"  of  the  surface,  and 
find  that  rates  of  crack  initiation,  early  growth  and  the  relaxation  of  re¬ 
sidual  surface  stresses  during  fatigue  are  all  sensitive  to  the  state  of 
ductility  of  the  surface.  The  depth  Into  the  surface  affected  by  environment 
is  too  shallow  to  permit  characterization  of  the  surface  ductility  by  conven¬ 
tional  techniques.  We  describe  the  development,  evaluation  and  refinement  of 
procedures  to  determine  the  ductility  of  the  surface  during  fatigue;  which 
rely  upon  characterization  of  the  microplasticity  of  the  surface,  in  plastic 
zones  at  the  tips  of  short  surface  cracks.  A  study  of  the  fundamental  mech¬ 
anisms  underlying  deformation  at  crack  tips  allows  us  to  compensate  for  the 
effect  which  crack  length  and  grain  size  has  on  the  microscopic  deformation, 
and  we  thus  obtain  fundamental  ductility  parameters  which  are  relatable  to 
alloy  type,  composition,  heat  treatment  and  humidity.  Five  parameters,  re¬ 
lated  to  different  aspects  of  surface  ductility,  are  defined  and  methods  to 
measure  them  are  discussed.  Results  of  an  initial  survey  of  the  relationship 
of  alloy  type  and  heat  treatment  to  the  sensitivity  of  surface  ductility  to 
humidity  are  also  described. 
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1.0  NOMENCLATURE 

Halfcrack  length 

Crack  tip  opening  displacement,  measured  at  a  tip  at  the  surface 
The  plastic  component  of  CTOO 
Grain  size 

Number  of  fatigue  cycles 

Number  of  fatigue  cycles  at  onset  of  Incubation 

Effective  maximum  surface  stress 

Distance  of  crack  tip  to  grain  boundary 

Material  parameter  related  to  critical  strain  energy  density 
required  for  propagation 

Opening  at  crack  center  at  stress,  crmax 

Effective  stress  Intensity  range 

Number  of  cycles  required  for  Incubation 

Strain  In  the  plastic  zone 

Normalized  strain  parameter  proportional  to 

Material  parameter  describing  rate  of  increase  In  strain  in  the 
plastic  zone 

Surface  shear  stress 

Surface  shear  stress  at  the  fatigue  limit 

Friction  shear  stress  on  cyclically  harden  surface 

Minimum  surface  shear  stress  producing  microplastic  deformation 

Maximum  tensile  stress  amplitude 

Alloy  yield  strength 

Coefficient  of  cyclic  hardening. 
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2.0  INTRODUCTION 

Environmental  humidity  substantially  affects  the  rate  at  which  cracks 
initiate  and  grow  on  the  surface  of  aluminum  alloys  during  fatigue.*’2  It 
appears  that  the  principal  effect  of  humidity  is  to  decrease  the  ductility  of 
the  alloy  surface.2  The  sensitivity  of  aluminum  alloys  to  humidity  varies 
with  alloy  type,  with  heat  treatment  and  apparently  with  minor  changes  in  the 
chemical  composition  of  an  alloy.  The  goal  of  this  research  is  to  model  the 
synergistic  effects  of  these  material  parameters  with  humidity  to  predict 
surface  ductility  and  hence  fatigue  behavior.  It  is  hoped  that  one  ultimate 
application  of  such  results  will  be  to  suggest  alloy  compositional  and  heat 
treatment  regimes  which  will  be  least  sensitive  to  fatigue  damage. 

Everyone  has  a  concept  of  the  meaning  of  ductility.  Ductility  as 
typically  encountered  in  metallurgy  is  not  a  fundamental  quantity  however,  but 
is  definable  only  in  terms  of  the  procedure  used  to  measure  it.  A  reduction 
in  area  test  simultaneously  assesses  the  contribution  from  several  factors 
which  affect  the  amount  of  deformation  a  material  will  sustain  before 
failure.  These  Include  hardenability,  and  resistance  to  crack  initiation  and 
to  crack  growth. 

The  purpose  of  the  first  year's  effort  in  this  program  was  to  develop 
and  evaluate  procedures  to  characterize  the  ductility  of  the  surface  of 
aluminum  alloys  subjected  to  fatigue.  Furthermore  the  plan  was  to  begin  a 
survey  of  the  effect  of  alloy  heat  treatment,  composition  and  humidity  on 
properties  of  ductility  using  the  new  experimental  techniques.  Our  intention 
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Is  to  define  a  set  of  measurements  which  characterize  surface  ductility  in  as 
fundamental  a  way  as  possible,  by  defining  procedures  to  measure  individual 
aspects  of  the  microplastic  deformation  properties  of  a  surface. 

•\ 

Early  growth  of  short  surface  cracks  involves  non-continuum  processes 
of  plastic  deformation  at  the  surface. ^  The  plastic  zone  size  of  short  cracks 
in  aluminum  alloys  is  essentially  determined  by  the  distance  of  the  crack  tip 
to  the  next  grain  boundary,  provided  that  the  applied  stress  is  kept  below  the 
alloy  yield  strength.  Slip  into  a  new  grain  only  begins  when  a  surface  crack 
tip  approaches  a  grain  boundary.  Our  measurements  of  "ductility"  rely  upon 
determination  of  the  rate  at  which  the  plastic  zone  develops  in  the  new  grain 
during  fatigue;  and  upon  determination  of  the  equilibrium  value  of  the  plastic 
component  of  opening  displacement  of  the  surface  crack  tip,  when  growth  into 
the  new  grain  begins  after  the  plastic  zone  is  mature. 

Interpretation  of  such  measurements  requires  a  fundamental  under¬ 
standing  of  microscopic  processes  involved  in  early  crack  growth.  Our  dis¬ 
coveries  In  the  past  year  regarding  growth  mechanisms,  their  relationship  to 
ductility  and  an  initial  survey  of  the  ductility  of  selected  alloys  are 
discussed  in  Section  4.0.  The  role  of  the  ductility  parameters,  we  define,  in 
crack  Initiation  and  in  relaxation  of  residual  surface  stresses  during  fatigue 
is  discussed  in  Section  5.0. 
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3.0  EXPERIMENTAL  PROCEDURES 

Early  crack  growth  measurements  are  made  using  miniature  flexural 
fatigue  specimens,  tapered  to  achieve  a  uniform  surface  stress  over  the 
specimen  gauge  section.  The  detailed  specimen  geometry  has  been  reported 
elsewhere.2  Specimens  are  prepared  by  machining  with  progressively  shallower 
cutting  depths  to  minimize  residual  surface  stress.  Final  cutting  depth  is 
approximately  20  urn  and  thereafter  specimens  are  mechanically  polished  ending 
with  a  0.05  urn  alumina  powder.  Residual  surface  stresses  obtained  with  this 
procedure,  and  measured  by  x-ray  diffraction,  are  small  (<20  MN/m2). 

Specimens  are  fatigued  in  laboratory  air  (humidity  measured)  and  are 
loaded  in  stroke  control  with  the  loading  fully  reversed  (R  =  0m-n/omax  =  -1) 
at  approximately  5  Hz.  They  are  then  transferred  at  intervals  to  either  a 
scanning  electron  or  optical  microscope  to  measure  the  crack  tip  deformation 
or  growth  parameters  of  interest. 

A  dry  box  has  been  constructed  which  can  be  used  to  maintain  the 
specimen  In  a  dry  air  environment.  This  is  flushed  for  several  hours  with  dry 
air,  and  then  a  cold  (liquid  nitrogen  cooled)  probe  is  exposed  to  pump  any 
residual  moisture.  We  have  found  the  system  workable  but  cumbersome  and  tine 
consuming  to  use.  To  replace  it,  we  have  evaluated  the  use  of  surface  coat¬ 
ings  which  could  prevent  penetration  of  ambient  moisture  to  the  surface  during 
fatigue.  Lacomit,  which  is  a  commercially  available  stop  coating  material 
used  in  electrooolishing,  was  found  to  provide  the  desired  characteristics. 

We  find  the  fatigue  crack  Initiation  and  early  growth  behavior  in  dry  air  and 
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with  specimens  coated  with  Lacomit  to  be  the  same.  By  coating  one  side  of  a 
specimen  with  Lacomit  and  leaving  the  other  exposed  to  ambient  humidity,  the 
effect  of  humidity  can  be  evaluated  on  a  single  specimen. 
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4.0  DUCTILITY  MEASUREMENTS 

Short  crack  growth  in  aluminum  alloys  is  quite  discontinuous.  Grain 
boundaries  can  substantially  impede  propagation  and  it  is  not  unusual  for 
short  cracks  to  be  completely  stopped  by  boundaries  over  the  entire  lifetime 
of  a  fatigue  specimen.  The  rate  at  which  cracking  proceeds  across  boundaries 
is  sensitive  to  the  ductility  of  the  surface,  and  so  the  measurement  of  appro¬ 
priate  parameters  associated  with  crack  growth  can  be  used  to  characterize 
ductility.  In  the  past  year  we  have  learned  that  it  requires  five  parameters 
to  specify  various  aspects  of  surface  ductility.  In  an  analogy  to  bulk 
measurements  these  can  be  thought  of  as:  (1)  the  cyclic  hardenability  of  the 
surface;  (2  and  3)  two  parameters  describing  the  micro-plasticity  of  the  sur¬ 
face  below  the  yield  strength,  a  magnitude  per  increment  strain  and  a  lower 
bound  friction  stress;  and  (4  and  5)  and  two  parameters  related  to  the  ulti¬ 
mate  material  strength.  It  appears  that  the  depth  below  the  surface  to  which 
the  ductility  parameters  can  be  altered  by  humidity  is  quite  shallow,  less 
than  a  few  micrometers.  Thus  conventional  procedures  can  not  be  used  to 
measure  the  parameters  of  Interest.  [An  attempt  was  made  to  use  Knoop  inden¬ 
tation  hardness  to  assess  the  change  of  surface  ductility  of  A1  2219-T851  with 
fatigue  In  humidity  vs  In  dry  air  with  negative  results].  While  the  data  are 
Incomplete,  It  also  appears  likely  that  three  of  the  parameters  are  sensitive 
to  humidity  as  well  as  to  the  material  variables.  Combinations  of  the  same 
parameters  are  required  to  describe  the  processes  of  crack  Initiation  and  the 
relaxation  of  residual  surface  stresses  In  aluminum  alloys,  as  well  as  the 
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growth  behavior  from  which  they  are  obtained.  Discussed  below  are  the  tech¬ 
niques  used  to  measure  the  ductility  parameters,  the  mechanisms  that  underlie 
the  behavior  characterized  and  the  status  of  our  survey  of  the  effect  of 
humidity  and  alloy  type  on  surface  ductility. 

4.1  History  of  Crack  Growth  Modeling  at  the  Science  Center 

We  have  been  involved  for  several  years  in  learning  how  to  measure 
the  closure  stresses  at  the  tips  of  small  surface  cracks.  This  ultimately  led 
to  the  development  of  an  empirical  procedure  to  estimate  the  closure  stress  at 
a  surface  crack  tip  by  measuring  the  opening  displacement  (CTOD)  at  the  max¬ 
imum  tensile  cyclic  stress  applied  during  fatigue.^  As  a  surface  crack  tip 
enters  a  new  grain,  a  closure  stress  develops  with  a  magnitude  proportional  to 
the  size  of  the  grain.  Closure  stresses  were  found  to  be  smaller  for  alloys 
fatigued  in  moist  than  in  dry  air,  and  were  found  to  vary  from  alloy  to  alloy 
for  a  fixed  grain  size.  It  was  shown,  by  comparing  predicted  and  measured 
rates  of  short  crack  growth  in  A1  2219-T851,  that  closure  stress  was  a  sub¬ 
stantial  factor  in  early  crack  growth.**5  The  initial  concept  upon  writing 
the  proposal  for  the  present  contract  was  to  use  CTOO  measurement  to  determine 
surface  ductility  variations  with  humidity  and  with  alloy  type.  Less  surface 
ductility  Is  associated  with  smaller  values  of  CTOO  and  hence  with  small 
closure  stresses. 

In  1979  a  series  of  experiments  were  done  under  the  support  of 
Rockwell  IR&D  funds  with  the  aim  of  predicting  early  growth  rates  from  the 
size  and  shape  of  grains  in  the  crack  path.  Grain  shape  parameters  were  used 
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to  calculate  the  closure  stress  at  a  crack  tip  at  the  surface,  giving  a 
predicted  rate  of  growth  based  upon  a  model. *  The  experiments  were  done 

using  A1  2219-T851  with  a  cyclic  stress  amplitude  of  +0.9  We  found 

that  variations  In  closure  stress  with  location  of  crack  tips  with  respect  to 
grain  boundaries  were  sufficient  to  describe  the  majority  of  the  growth  rate 
observations,  with  one  exception.  Propagation  rates  slow  down  just  before 
growth  begins  Into  a  new  grain,  as  slip  commences  across  the  boundary.  We 
found  that  the  retardation  In  growth  rate  was  larger  than  predicted  from  the 
closure  model  If  a  crack  tip  encountered  a  small  grain,  and  that  with 
sufficiently  small  fatigue  increments  one  could  see  that  propagation  actually 
stopped.6 

In  1978  Chang  _et  jil,7  had  proposed  a  model  of  early  crack  growth  which 
considered  the  dynamics  of  development,  during  fatigue,  of  a  plastic  zone  at  a 
crack  tip  In  an  initially  undeformed  medium.  It  was  hypothesized  that  de¬ 
velopment  of  a  mature  plastic  zone  required  accumulated  microplastic 
deformation,  and  it  was  shown  that  In  this  case  crack  growth  would  be  arrested 
until  a  critical  strain  energy  density  was  developed  at  a  crack  tip.  The 
theory  was  first  used  to  model  the  duration  of  non-propagation  from  Inter- 
metallic  particles  fractured  during  fatigue  of  A1  2219-T851,  and  successfully 
predicted  observed  statistical  variations  In  non-growth  periods  with  particle 
size  and  with  encompassing  grain  size.8 

It  was  realized  that  the  Chang  model  had  the  type  of  dependence  upon 
slip  distance  required  to  explain  the  Inaccuracy  of  our  short  crack  growth 
model.  The  applicability  of  the  Chang  model  to  crack  growth  Into  a  new  grain 

8 

C2661A/bw 


i 


Rockwell  International 

Scwrtc*  Center 

SC5211.1TR 


was  easily  rationalized,  as  we  knew  that  slip  into  the  new  grain  did  not  occur 
until  just  before  a  crack  tip  entered  a  new  grain. ^  Thus  growth  across  a 
boundary  begins  essentially  into  an  undeformed  medium.  We  refer  to  the  pro¬ 
cess  during  which  the  plastic  zone  develops  as  incubation.  It  was  found  that 
a  refurbished  model,  incorporating  the  effects  of  grain  size  and  shape  on  both 
closure  stress  and  incubation  period,  was  quite  adequate  to  predict  our  ex¬ 
perimental  growth  observations.4’® 

At  this  point  it  was  quite  clear  that  the  static  model  of  CTOD  be¬ 
havior  upon  which  our  closure  stress  predictions  were  based  was  too  sim¬ 
plistic.  Fortunately,  owing  to  its  strong  dependence  on  cyclic  stress 
amplitude,  incubation  was  only  of  secondary  importance  to  crack  growth  for  the 
large  amplitudes  used  in  the  growth  experiments  (±0.9  oyield)  allowing  the 
effect  of  closure  stress  to  be  independently  assessed.  We  expected  incubation 
could  be  the  dominant  factor  in  controlling  retardation  of  short  crack  growth 
across  grain  boundaries  for  low  cyclic  stress  amplitude,  and  that  the 
extension  of  the  duration  of  the  Incubation  period  at  low  stress  amplitude 
would  allow  us  to  assess  the  dynamics  of  the  microplasticity  at  a  surface 
crack  tip. 

At  about  the  same  time  a  new  aspect  of  the  effect  of  humidity  on 
surface  ductility  was  discovered.  An  ongoing  project  supported  by  Rockwell 
has  been  directed  at  refining  and  testing  a  model  (Initially  developed  under 
ONR  Contract  No.  00014-76-C-0452)  of  fatigue  crack  initiation  at  constituent 
surface  particles.  We  found  that  particles  were  more  easily  fractured  by 
fatigue  in  dry  than  in  moist  air,  and  furthermore  that  our  simple  particle 
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fracture  model  substantially  overestimated  the  numbers  of  particles  actually 
broken  at  low  cyclic  stress  amplitudes,  particularly  later  in  the  fatigue 
lifetime.  A  new  initiation  model  was  proposed  in  which  we  assumed  that  the 
friction  stress  (the  minimum  shear  stress  required  for  microplastic  deforma¬ 
tion)  Increased  progressively  during  fatigue,  at  a  rate  sensitive  to  the 
environmental  humidity.  A  simple  model  of  the  process  was  found  adequate  to 
explain  our  experimental  observations.0 

It  was  recognized  that  such  cyclic  hardening  of  the  surface  must  also 
effect  the  Incubation  mechanism  by  altering  the  rate  of  slip  Into  a  new  grain. 
The  Incubation,  cyclic  hardening  and  closure  stress  phenomena  are  sensitive  to 
different  aspects  of  the  ductility  of  an  alloy  surface.  The  major  thrust  of 
our  research  on  this  contract  In  the  past  year  has  been  directed  at  modeling 
these  processes,  and  in  definition  of  procedures  to  extract  meaningful  ductil¬ 
ity  parameters  from  such  measurements.  The  major  results  of  our  research  In 
the  past  year  are  as  follows. 

1.  The  physical  reality  of  Incubation  has  been  demonstrated  by 
measuring  the  progressive  Increase  In  the  plastic  component  of 
CTOD  (CTOOp)  for  crack  tips  arrested  at  grain  boundaries. 

2.  We  find  that  cyclic  hardening  of  the  surface  Increases  the 
duration  of  Incubation  by  suppressing  slip  into  a  new  grain. 
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3.  A  model  of  the  incubation  period  has  been  constructed  which 
considers  the  effect  of  grain  size,  crack  length,  cyclic 
hardening,  grain  crystallographic  orientation  and  stress 
amplitude  on  the  incubation  period.  Predictions  of  several 
aspects  of  this  model  have  been  tested  and  the  model  appears  to 
be  good. 

4.  A  model  of  the  effect  of  constraint  of  slip  by  grain  boundaries 
on  closure  stress  has  been  developed,  whose  predictions  agree 
with  our  empirical  closure  stress  equations.  A  simplified  model 
of  the  progressive  increase  of  CT0Dp  during  incubation  has  also 
been  devised  and  appears  adequate  to  explain  our  observations. 

5.  A  set  of  five  parameters  describing  surface  ductility  and 
procedures  to  measure  these  parameters  have  been  defined. 

6.  An  initial  survey  of  the  effect  of  heat  treatment  and  humidity 
on  surface  ductility  is  complete. 

4.2  Equilibrium  Growth  Phenomena 

Rice9  has  suggested  that  the  plastic  component  of  crack  tip  opening 
displacement  CT0Dp  -  ey z0/2,  where  zp  Is  the  size  of  the  plastic  zone  and  ey 
is  the  strain  at  the  yield  strength.  For  short  surface  cracks  In  aluminum 
alloys  fatigued  below  the  yield  strength,  zQ  is  the  distance  of  the  cra„k  tip 
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to  the  next  grain  boundary.-*  A  measurement  of  CT0Dp  vs  zQ  can  therefore  be 

used  to  obtain  ey,  a  parameter  which  describes  one  facet  of  ductility. 

Results  typical  for  A1  2219-T851  alloy  are  shown  in  Fig.  1.  Parameters  of 

short  surface  cracks  are  measured  in  the  SEM  using  a  jig  which  permits 

specimens  to  be  loaded  to  the  maximum  tensile  stress  applied  during  fatigue 

(Fig.  2).  Measurements  made  for  each  crack  tip  are  illustrated  in  Fig.  3.  We 

infer  from  results  of  an  analysis  by  Green  and  Sneedon1^  that  for  an  elastic 

crack  (i.e.,  as  zQ  +  0)  CT0D/6(omax)  -  0.1.  The  parameter  6 ( ama x }  is  the 

opening  at  the  center  of  the  crack  and  is  approximately  proportional  to  the 

crack  length  2c.  Thus  with  the  normalization  used  in  Fig.  l,  all  deformation 

above  the  dashed  line  is  plastic.  For  convenience,  we  define  a  normalized 

strain  given  by  Eq.  (D,  which  represents  a  normalized  failure  strain  for 

the  material  in  environment.  The  actual  strain  e  =  0.0182  o  e'/o  .  . 

y  max  y  yield 

The  proportional ty  coefficient  comes  from  the  ratio  of  crack  opening  to  crack 
length,  which  is  only  approximately  known.  Thus  e y  is  determinable  with  more 
accuracy  than  is  ey. 


2c 

z_ 


CTOD 

ZTo - 

max 


-  0.1 


(1) 


Values  of  ey  measured  in  this  manner  for  alloys  A1  2219-T851  and  A1 
7075-T6  are  listed  in  Appendix  A.  Attempts  to  measure  for  A1  2024-T3 
yielded  the  erroneous  result,  ey  *  0.  This  is  an  artifact  of  the  measurement 
procedure  which  arises  when  the  cracks  in  the  data  set  are  crystallographic. 
Apparently  the  plastic  deformation  associated  with  the  growth  of 


is 
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Fig.  1  Crack  tip  opening  displacement  vs  distance  z0  of  the  crack 
tip  to  the  next  grain  boundary  for  A1  2219-T851.  With  the 
normalization  used  a  value  of  CT00/6(amax)  *  0.1  is  the 
elastic  component  of  the  opening  displacement.  Plastic 
component  Increases  with  increased  distance  of  crack  tip 
to  boundary.  Dashed  line  Is  predicted  elastic  opening 
displacement. 
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Dimensional  parameters  used  in  determination  of  microplastic 
deformation.  CTOO  is  found  by  extrapolation  of  displacements 
made  near  the  tip  to  the  tip  as  explained  elsewhere.3 
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crystallographic  cracks  Is  entirely  In  the  plane  of  the  crack.  A  CTOD 
measurement  on  such  cracks  detects  only  the  elastic  opening  normal  to  the 
crack  path  with  the  result  that  CT00/6(omax)  -  0.1,  independent  of  z0«  Cracks 
which  obviously  were  crystallographic  were  excluded  from  the  data  set  in  the 
2219  and  7075  alloys.  An  extensive  survey  of  surface  cracks  in  A1  2024-T3  did 
not  locate  any  non-cry stall ographic  cracks.  A  possible  alternative  for 
crystallographic  cracks  is  to  measure  displacements  in  the  plane  of  the 
crack.  Such  an  analysis  is  in  progress. 

As  an  interim  step  in  development  of  a  model  of  the  dynamics  of 
microplastic  deformation  at  a  surface  crack  tip,  we  assessed  approaches  to 
predict  the  equilibrium  CT0Dp  and  closure  stress  values  encountered  when 
plastic  zone  size  is  constrained  by  the  grain  boundaries.  In  Fig.  4  closure 
stress  values  predicted  as  a  function  of  distance  of  crack  tip  to  grain 
boundary,  using  two  different  models,  are  compared  to  experimental  data.  One 
prediction  is  obtained  using  a  finite  element  analysis  of  the  deformation  at  a 
crack  tip  with  the  plastic  zone  size  constrained  by  a  boundary.  The  second 
approach  involves  use  of  a  superdislocation  model  of  the  deformation 
process.  Details  on  both  models  are  provided  in  Appendix  B.  A  crack  closure 
stress  monotonlcally  increasing  with  increased  distance  of  the  crack  tip  to 
the  grain  boundary  is  both  observed  and  predicted.  It  appears  that  there  is  a 
tendency  for  some  of  the  closure  stress  values  to  occur  at  systematically 
smaller  values  than  expected  with  either  model.  We  believe  this  is  a  result 
of  the  dynamics  of  plastic  zone  maturation  as  discussed  next. 
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Fig.  4  Ratio  of  closure  stress  to  maximum  tensile  stress  increases 
with  increased  distance  U)  of  crack  tip  to  the  next  grain 
boundary.  Here  (a)  Is  crack  depth  in  a  two  dimensional 
crack  model . 
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4.3  Dynamics  of  Plastic  Zone  Maturation 

Once  slip  from  a  surface  crack  tip  begins  into  a  new  grain,  crack 
growth  stops. The  duration  of  the  incubation  period  is  longer  the  smaller 
the  cyclic  stress  amplitude.  Using  A1  2219-T851,  a  series  of  experiments  were 
done  to  assess  the  dynamics  of  the  incubation  process  by  measuring  CTOD  vs 
fatigue  cycles  N  for  a  large  number  of  crack  tips.  Specimens  were  prefatigued 
at  omax  =  0.9  Oyieu  to  generate  a  distribution  of  short  cracks.  The  cyclic 
stress  amplitude  was  then  lowered  and  the  crack  growth  monitored  at  increments 
during  fatigue.  In  particular  we  were  hoping  to  locate  cracks  which  initially 
continued  to  grow,  but  soon  encountered  grain  boundaries  at  the  lower  stress 
amplitude  and  stopped.  After  an  experiment  was  complete  the  surface  was 
chemically  etched  to  reveal  the  location  and  size  of  grains  in  the  path  of 
each  crack  tracked. 

In  Fig.  5  we  show  three  results  typical  of  such  experiments.  The 
data  are  plotted  using  the  same  format  as  in  Fig.  1,  and  show  a  progressive 
increase  in  CTOD  during  fatigue,  which  illustrates  the  dynamics  of  accumula¬ 
tion  of  plastic  deformation  as  a  crack  tip  transits  a  grain  boundary.  In  case 
A,  both  crack  tips  are  undergoing  incubation,  and  the  crack  length  (2c)  is 
constant  during  fatigue.  The  progressive  increase  In  CTOD  during  fatigue 
parallels  an  increasing  density  of  dislocations  in  the  next  grain  with  each 
fatigue  cycle.  Case  B  differs  from  case  A  only  in  that  the  other  surface 
crack  tip  was  not  at  a  boundary  and  was  propagating  so  that  2c  increased 
progressively  during  fatigue.  In  case  C,  the  applied  stress  was  less  than  the 
friction  stress  required  for  deformation  to  begin  Into  the  next  grain  and  the 
CTOD  remained  elastic.  [We  find  that  in  the  latter  case  ultimate 
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Fig.  5  Three  examples  of  the  progressive  change  in  CTOO  during 
incubation,  for  c^cV  tips  at  grain  boundaries. 
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propagation  is  typically  intergranular  and  results  from  the  fracture  of  grain 
boundaries  ahead  of  but  not  directly  connected  to  the  main  crack]. 

We  propose  expressions  in  Section  4.3.1  and  4.3.2  to  predict: 


1.  The  duration  of  the  incubation  period,  AMj,  from  which  a 

coefficient  0  describing  the  hardenabil ity  of  the  surface  is 
determined. 


2.  The  increase  in  the  plastic  component  (CT0Dp)  of  CTOD  during 
fatigue. 


4.3.1  The  Duration  of  the  Incubation  Period  ( ANT )  Before  Commencement  of 


/ANl+Ns  5  D/ScTnT  (t  -  x  )2e"26NdN  =  6.  (2) 

..  eft  o 

NS 

Eq.  (2)  is  obtained  from  an  analysis  of  the  incubation  process  by  Chang  et 
al ,7  with  several  empirical  modifications.  S  and  6  are  material  parameters. 
These  are  related  to  alloy  ductility,  as  discussed  later.  D  is  the  slip 
distance  into  the  new  grain  and  2c (N)  is  the  crack  length  after  N  fatigue 
cycles  (Fig.  3).  N$  is  the  number  of  fatigue  cycles  at  which  incubation 
began.  We  define  an  effective  surface  shear  stress  (Teff  -  T0)e'eN  to  replace 
the  local  stress  S  In  the  Chang  derivation.  Due  to  cyclic  hardening  of  the 
surface  the  friction  stress  (x^ )  at  the  surface  is  assumed  to  increase  at  a 
rate  proportional  to  the  microplastic  deformation  at  the  surface.  Thus 
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eff 
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(3) 


S  =  (Teff  -  V  s  (Teff  “  To)e 


-0N 


(4) 


0  is  a  coefficient  of  hardenability.  rQ  is  the  friction  stress  in  unfatigued 
material.  is  an  effective  shear  stress  which  is  determined  by  the 

crystallographic  orientation  of  the  grain  at  the  crack  tip.  re^  is  con¬ 
structed  for  a  given  grain  orientation  by  calculating,  from  the  applied 
stress,  the  shear  stress  in  each  of  the  (111)  slip  planes  in  the  fee  alumi¬ 
num.  The  largest  component  of  these,  lying  both  in  a  (111)  plane  and  in  the 
surface,  is  defined  to  be  (See  Figs.  6  and  7). 

In  the  special  case  that  the  crack  length  is  constant  during 
incubation  and,  if  0  is  small,  Eq.  (2)  can  be  rewritten  as 

InUNjD/Tc}  +  ln[U/0)(Teff  -  tq)2]  =  20N$  .  (5) 

A  series  of  experiments  were  done  to  test  Eq.  (5),  from  which  we 
could  determine  0  for  selected  alloys.  Specimens  were  prefatigued  in  dry  air 
to  produce  a  set  of  surface  cracks  with  a  maximum  length  of  about  150  urn. 
Several  micrometers  of  the  surface  were  then  removed  by  mechanical  polishing, 
to  return  the  surface  to  an  unhardened  state.  The  surface  was  then  chemically 
etched  to  reveal  the  grain  boundaries,  and  duration  of  incubation  periods  for 
crack  growth  across  grain  boundaries  was  determined  during  subsequent 
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fatigue.  This  was  accomplished  by  making  replicas  of  the  surface  at  small  in¬ 
tervals  during  fatigue,  which  were  later  examined  to  measure  N^,  D,  2c  and 
ANj.  Results  for  A1  2024-T3  for  two  humidities  are  shown  in  Fig.  8.  0  is 

found  from  the  slope  of  the  incubation  plot.  Results  for  all  materials 
studied  to  date  are  given  in  Appendix  A.  The  effect  of  humidity  is  to  in¬ 
crease  the  magnitude  of  0.  TEM  of  the  dislocation  structure  near  the  surface 
of  A1  2024-T3  show  dislocation  tangles  in  the  alloy  exposed  to  humid  air, 
while  zones  of  persistent  slip  are  found  on  a  specimen  fatigued  in  dry  air 
(Fig.  9). 


4.3.2  CT0Dr  vs  N 

We  obtained  the  theoretical  dependence  of  CT0Dp  on  number  of  fatigue 
cycles  N  using  the  model  reported  in  Ref.  7,  and  by  assuming  that  CT00p  is 
proportional  to  the  number  of  dislocations  In  the  slip  band  at  the  crack 
tip.  Thus 


CT0Dp  *  02/Ic  (H  -  Ns)(xeff  -  io)2e_2eN  .  (6) 

We  have  assumed  that  the  crack  length  Is  constant  during  incubation.  An 
experimental  test  of  Eq.  (6)  using  A1  2219-T851  is  shown  In  Fig.  10.  The 
progressive  Increase  In  CT0Dp  with  fatigue  Illustrates  the  physical  reality  of 
the  Incubation  process.  Numerous  examples  of  the  effect  have  been  studied  of 
which  Fig.  10  is  typical. 
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Fig.  10  Progressive  Increase  in  CT0DQ  for  crack  tip  which 
reached  a  boundary  at  15  *  103  fatigue  cycles. 
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5.0  MEW  INSIGHT  INTO  SURFACE  DUCTILITY 

The  ductility  of  the  surface  of  aluminum  alloys  changes  progressively 
during  fatigue.  Surface  hardening  occurs  both  in  dry  and  in  humid  air,  but  is 
much  more  rapid  and  greater  In  magnitude  in  humid  air.6  The  effect  of  the 
hardening  is  seen  both  in  the  crack  initiation  and  in  the  early  growth  behav¬ 
ior  of  the  alloys.  For  instance,  in  A1  2219-T851  humidity  decreases  the  num¬ 
bers  of  intermetallic  particles  which  fracture  during  fatigue.  Furthermore, 
humidity  decreases  the  rates  at  which  residual  surface  stress,  induced  by 
machining  in  A1  2219-T851,  relax  during  fatigue. ^  [Figures  11  and  12  which 
illustrate  the  effect  of  cyclic  hardening  of  the  surface  on  particle  fracture 
and  relaxation  of  residual  stress  were  obtained  in  a  Rockwell  IR&D  supported 
investigation.]  Both  these  effects  are  attributed  to  the  reduction  in 
microplastic  deformation  per  cycle  which  results  from  an  environmentally 
assisted  increase  in  friction  stress  of  the  near  surface  material. 

During  early  crack  propagation  the  effect  of  ductility  on  growth  rate 
is  complex.  Slip  which  promotes  propagation  for  fatigue  amplitudes  below  the 
yield  strength  does  not  extend  In  substantial  amounts  beyond  a  grain  boundary 
until  a  surface  crack  tip  approaches  the  boundary.  Propagation  ceases  at  the 
boundary  until  a  mature  plastic  zone  Is  developed  in  the  new  grain.  We  refer 
to  this  process  as  incubation.  Propagation  rate  after  incubation  is  complete 
Is  controlled  by  the  closure  stress  at  the  surface  crack  tip,  which  is  a 
function  of  the  grain  size  and  of  the  ductility  of  the  alloy  at  the  surface. 
The  effect  of  a  decreased  ductility  during  fatigue  in  humid  air  is  to  increase 
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Fig.  11  Density  of  particles  fractured  on  the  surface  of  A1  2219-T851 
during  fully  reversed  fatigue  for  four  maximum  cyclic  stress 
amplitudes.  Symbols  are  experimental  data  obtained  by  SEM. 
Curves  are  computer  predictions.  Decrease  rate  of  fracture 
later  in  lifetime  is  due  to  cyclic  hardening  of  the  surface 
and  is  more  pronounced  for  fatigue  in  humid  than  in  dry  air. 
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the  Incubation  period  and  to  decrease  the  closure  stress,  leading  to  accel¬ 
erated  rates  of  propagation  after  incubation.  The  net  effect  of  ductility  on 
lifetime  is  therefore  determined  by  the  distribution  in  size  of  the  grains  at 
the  alloy  surface. 

We  also  find  that  surface  hardening  has  no  appreciable  effect  on 
crack  nucleation  at  grain  boundaries.  We  conclude  that  slip  leading  to  crack¬ 
ing  of  particles  and  slip  which  accompanies  microcrack  propagation  is  long 
range,  extending  across  the  grain  matrix  to  the  grain  boundaries.  Grain 
boundary  cracking  involves  void  nucleation,  requiring  only  short  range 
deformation.12  Models  based  upon  this  view  have  been  highly  successful  in 
predicting  crack  initiation  and  growth  behavior  in  A1  2219-T851.4*6’8  We 
believe  that  this  situation  pertains  to  high  strength  aluminum  alloys  in 
general . 

The  effect  of  humidity  is  to  reduce  the  propensity  for  long  range 
slip.  TEM  of  near  surface  dislocation  structure  shows  banded  slip  in  an 
A1  2024-T3  alloy  fatigued  in  dry  air  and  dislocation  tangles  in  humid  air 
(Fig.  9).  Surface  hardening  is  treated  theoretically  by  the  introduction  of  a 
friction  shear  stress  tq,  altered  during  fatigue  at  a  rate  determined  by,  9,  a 
hardenability  coefficient. 

Our  approach  to  measuring  surface  ductility  has  been  to  recognize 
that  early  crack  growth  behavior  is  affected  both  by  geometrical  factors 
•  (grain  size,  crack  length,  crack  location)  and  by  the  ductility  of  the 
surface.  We  have  formulated  relationships  through  combinations  of  models  and 


31 

C2661A/bw 


Rockwell  International 

Science  Center 
SC5211.1TR 


empirical  observations  which  allow  us  to  relate  measurables  to  the  geometrical 
parameters  and  to  the  desired  ductility  parameters.  By  this  means  we  find 
that  five  distinct  material  parameters  are  required  to  describe  the  behavior 
of  an  aluminum  alloy  surface,  as  affected  by  humidity.  These  are: 

1.  Parameter  0,  which  describes  the  increase  in  friction  shear 
stress  t0,  Airing  fatigue. 

2.  Parameter  £,  which  specifies  the  rate  of  increase  in  dislocation 
density  in  a  slip  band  at  a  crack  tip  with  cyclic  strain. 

3.  Parameters  6  and  e^,  which  are  related  to  a  critical  strain 
energy  density  for  propagation  and. an  ultimate  material  strength 
respectively. 

All  five  parameters  are  independent  of  N  and  of  alloy  grain  size,  and 
two  are  known  to  be  functions  of  humidity.  (By  defi nation  tq  is  independent 
of  humidity,  since  any  humidity  dependence  is  carried  in  0).  0  is  larger  and 

is  smaller  in  moist  air  than  in  dry  air.  The  effect  of  humidity  on  e  and  5 
remains  to  be  established,  as  do  the  effects  of  alloy  composition  and  heat 
treatment  on  these  parameters.  Initial  results  suggest,  however,  that  0  is 
much  more  sensitive  to  humidity  in  underaged  alloys. 

The  incubation  model  is  based  upon  a  "small  deformation"  theory.  The 
Increase  in  the  friction  stress  during  fatigue  has  been  incorporated 
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empirically  into  the  model  to  explain  the  progressive  hardening  of  the  surface 
manifested  in  both  crack  initiation  in  early  growth  behavior  observed  in  2000, 
6000  and  7000  series  alloys.  Our  superdislocation  model  of  CT0D  relies  upon  a 
large  deformation  analysis  to  predict  equilibrium  values  of  CTOD  and  of  crack 
tip  closure  stress.  These  two  methods  both  predict  that  a  critical  strain 
energy  density  at  the  crack  tip  is  required  for  propagation.  We  denote  this 
effect  by  the  parameter  0  for  the  incubation  model  a the  parameter  for 
the  CTOD  model.  The  interrelationship  between  0  and  remains  to  be  deter¬ 
mined.  The  cycle  dependence  of  CTOD  during  incubation  is  also  potentially 
rich  in  information  about  surface  ductility,  but  our  new  model  of  this  process 
must  be  tested  before  it  can  be  generally  applied.  The  status  of  what  can 
presently  be  measured  regarding  ductility  is  summarized  in  Table  1. 
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Table  1 


Summary  of  Ductility  Measurements 


Measurement 

Process/ 

Parameter 

Incubation 

Period 

CTOD  vs  N 
Dependence 

CTOD 

Equilibrium  Value 

Increase  in 

T0/(e) 

Measured  using 
Eq.  (5) 

Measured  using 
Eq.  (6) 

Not  expected  to 
directly  affect 
equilibrium  CTOD 

Rate  of  increase 
of  dislocation 
density  with 
cyclic  strain/ 

(C) 

The  ratio  C/0 
measured  using 
Eq.  (5) 

e 1 

measured 

P 

using  Eq.  (6) 

Not  expected  to 
directly  affect 
equilibrium  CTOD 

Critical  strain 
energy  density 
for  propagation/ 

(0,  ey) 

The  ratio  C/0 
measured  using 
Eq.  (5) 

e' 

— measured 

P 

using  Eq.  (6) 

e'  measured  using 
Eq*  (1) 

The  three  material  parameters  C,  (0,  and  e  each  describe  a 
different  facet  of  ductility,  but  it  is  likely  that  they  are  partially 
interrelated.  This  prospect  will  require  a  further  improvement  in  understand¬ 
ing  of  the  nature  of  near  surface  deformation  in  aluminum  alloys.  Hopefully, 
in  the  second  year  of  research,  a  unified  theory  of  surface  ductility  will  be 
forthcoming.  Intuitively,  each  parameter  represents  a  different  aspect  of 
ductility  such  as  might  be  encountered  in  a  macroscopic  evaluation  of  alloy 
ductility.  (3,  Cy)  is  similar  to  a  large  deformation,  reduction  in  area 
parameter.  Materials  with  large  (0,  e^)  are  more  ductile  and  will  sustain 
more  strain  at  a  microcrack  tip  before  separation.  Hardening  or  softening  of 
the  bulk  material  during  fatigue  is  common  and  can  arise  from  dislocation 
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mechanisms.  Parameter  6  describes  the  rate  of  surface  hardening,  which  is 
highly  sensitive  to  presence  of  humidity  for  aluminum  alloys.  The  depth  of 
material  at  the  surface  affected  by  environment  is  very  shallow,  on  the  order 
of  several  microns,  as  suggested  by  the  amount  of  material  which  must  be 
removed  to  return  the  surface  to  an  unhardened  state.  Parameters  £  and  tq  are 
distantly  related  to  the  microplastic  stress-strain  response  of  a  material, 
which  can  be  determined  for  bulk  material.  Because  the  slip  important  to 
early  crack  growth  is  confined  to  a  thin  layer  near  the  surface,  however, 

(B,  £y),  0  and  £  are  all  sensitive  to  humidity  and  to  the  partial  relaxation 
of  bulk  constraints  on  deformation  at  the  surface. 
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6.0  DISCUSSION  ANO  SUMMARY 

It  has  been  many  years  since  Frost*"*  first  concluded  that  the  plain 
fatigue  limit  was  a  function  of  grain  size  and  that  crack  propagation  across 
grain  boundaries  required  lower  cyclic  stress  amplitudes  the  larger  the  grain 
hence  the  larger  the  initial  crack  size.  Forrest  and  Tate*4  determined  that 
an  empirical  expression  proposed  by  Frost,  which  related  crack  length  to  the 
cyclic  stress  amplitude  required  for  propagation,  required  only  a  slight 
modification  to  describe  the  effect  of  grain  size  on  the  fatigue  limit  of  a 
70/30  brass.  Later  Hoeppner*5  showed  that  for  stress  amplitudes  above  the 
fatigue  limit  a  smaller  grain  size  suppressed  propagation  across  grain 
boundaries  for  more  fatigue  cycles  than  a  large  grain  size.  Thus  Hoeppner 
observed  the  process  we  refer  to  as  incubation.  Taira,  Tanaka  and  Nakai*6 
have  proposed  a  model  of  the  fatigue  limit  based  upon  a  derivation  which  shows 
that,  below  a  critical  stress  amplitude,  slip  from  a  crack  tip  does  not 
proceed  across  a  boundary  into  the  next  grain.  By  assuming  propagation  across 
the  boundary  requires  slip  into  the  new  grain,  they  derive  a  dependence  of  the 
fatigue  limit  upon  grain  size  which  appears  adequate  to  explain  experimental 
observations  of  the  effect  of  grain  size  on  fatigue  limit  in  a  low  carbon 
steel. *7  The  model  of  the  incubation  process  proposed  by  Chang  et  al.^ 
distinguishes  the  effects  of  crack  length  from  that  of  grain  size  at  the  crack 
tip.  It  differs  from  the  model  In  Ref  16  In  that  it  addresses  crack 
propagation  as  well  as  the  fatigue  limit.  It  presumes  that  slip  into  the  next 
grain  is  necessary  but  not  sufficient  for  propagation,  and  that  in  addition 
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sufficient  microplastic  deformation  must  accrue  for  propagation.  The  Chang 
model  predicts  a  fatigue  limit  (t^)  for  intergranular  crack  propagation  of  the 
form 

t e  -  constant  /07TT  +  t  .  (7) 

f  o 

There  are,  thus,  three  separate  equations  [Eq.  (2),  (5)  and  (7)] 
related  to  various  aspects  of  the  incubation  process,  which  require 
evaluation.  Of  these  Eq.  (2),  and  its  approximate  representation  Eq.  (5),  has 
been  subjected  to  the  most  scrutiny.  The  effect  of  grain  size,  crack  length 
and  hardenabillty  on  the  duration  of  the  incubation  period  appears  to  be 
correct.  The  dependence  upon  stress  amplitude  is  known  to  be  approximately 
correct  but  has  not  been  tested  rigorously.  The  effect  of  crystallographic 
orientation  on  the  Incubation  period  has  been  verified  using  a  statistical 
analysis  but  not  by  direct  measurement  of  the  crystallographic  orientation  of 
a  grain  at  a  crack  tip.®  We  know  the  duration  of  incubation  is  correctly 
described  by  Eq.  (2)  if  propagation  Into  the  new  grain  is  transgranular  and 
either  crystallographic  on  non-crystallographic.  We  do  not  know,  however,  if 
(0/C)  is  the  same  for  both  crystallographic  and  noncrystallographic  cracks  in 
the  same  alloy.  Furthermore  we  have  not  looked  for  a  possible  effect  of  crack 
depth  (as  compared  to  crack  length)  on  the  duration  of  Incubation. 

In  regard  to  prediction  of  CT00p  [Eq.  (6)]  only  the  dependences  upon 
0  and  N  have  been  tested  and  these  appear  correct.  It  can  be  argued  that, 
since  Eqs.  (2)  and  (6)  come  from  the  same  root  equation,  their  dependencies 
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upon  D  and  2c  should  be  equally  valid.  However,  we  would  not  be  surprised  if 
this  were  not  the  case,  considering  the  many  assumptions  and  approximations  in 
the  theory.  Evaluation  of  Eq.  (6)  is  extremely  tedious  owing  to  the 
difficulty  of  measuring  CTOOp  vs  N.  In  the  alloys  studied  to  date,  Eq.  (6)  is 
invalid  for  crack  lengths  longer  than  *  300  urn,  for  which  length  deformation 
at  the  crack  tip  is  determined  by  continuum  constraints. 

In  regard  to  Eq.  (7),  we  have  serious  doubts  of  its  applicability  to 
aluminum  alloys.  The  2000  series  aluminums  on  which  we  have  consentrated  the 
majority  of  our  research  do  not  have  an  endurance  limit.  At  low  stress 
amplitudes  crack  growth  at  the  surface  is  intergranular  and  Eq.  (7)  can  not  be 
expected  to  apply.  The  growth  process  in  this  case  involves  fracture  of  grain 
boundaries  ahead  of  but  not  directly  connected  to  the  crack  tip,  followed  by 
linking  of  the  secondary  cracks  with  the  main  crack.  We  have  measured  the 
CTOD  of  the  main  crack  and  find  that  intergranular  cracking  can  occur  even  if 
the  cyclic  hardening  has  proceeded  to  the  point  that  the  CTOD  remains 
elastic.  It,  thus,  appears  that  the  fracture  of  grain  boundaries  ahead  of  the 
crack  tip  can  result  from  strain  amplification  due  to  the  elastic  strain 
associated  with  the  main  crack.  Perhaps  Eq.  (7)  might  apply  to  alloys  which 
have  a  true  endurance  limit  such  as  certain  steels. 

In  summary,  the  models  from  which  ductility  parameters  have  been  de¬ 
fined  and  calculated  are  not  completely  tested.  Indeed,  given  the  multi¬ 
plicity  of  micros true tural  parameters  Involved  it  is  not  likely  that  at  a 
complete  evaluation  of  the  equations  will  be  forth  coming  under  this  con¬ 
tract.  Furthermore,  the  mechanism  of  hardening  of  the  surface  is  not 
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presently  known,  although  we  suspect  that  It  may  result  from  the  effect  of  the 
oxide  layer  on  the  dislocation  structure  near  the  surface.  Grosskreutz18  has 
reached  such  a  conclusion  for  A1  1100  alloy.  However,  the  ductility 
parameters  have  demonstrated  utility  in  predicting  several  aspects  of  fatigue 
behavior.  The  parameters  must  be  viewed  as  empirical  quantities  until  model 
verification  is  complete  and  mechanisms  are  further  defined.  Until  then  they 
serve  as  useful  evaluators  of  the  sensitivity  of  materials  to  fatigue  failure, 
and  guide  in  determination  of  the  effect  of  alloy  chemistry  and  heat  treatment 
on  fatigue  behavior  in  humid  air. 
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APPENDIX  A 

EXPERIMENTAL  DUCTILITY  PARAMETER  VALUES 


Alloy 

Humidity,  % 

ey 

0,  cycles"1 

e/C,MN//m 

Tp.MN/m2 

A1 

2215-T851 

0 

0.36 

1  x  io-5* 

- 

182 

A1 

2215-T851 

60 

0.25 

3  x  io-5 

CM 

1 

O 

r-i 

X 

e 

iH 

182 

A1 

2024-T3 

0 

- 

4.2  x  io-5 

0.9  x  IO"2 

- 

A1 

2024-T3 

27 

- 

3.1  x  io-4 

0.9  x  IO"2 

- 

A1 

2024-T6t 

50 

- 

<1  x  10"5 

2  x  io-1 

- 

A1 

7075-T6 

50 

0.40 

3.5  x  IO"5 

rH 

1 

o 

rH 

X 

r^ 

*H 

- 

*A  second  heat  of  AT  2219-T851  has  been  found  for  which  9  <<  l  x  io-5 
and  is  independent  of  humidity.  The  relationship  of  e  to  alloy  com¬ 
position  is  now  being  assessed. 
tCommercial  T3  material,  aged  to  peak  hardness. 
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ON  CRACK  CLOSURE  PHENOMENON  IN  FATIGUE 


Roger  Chang 

Rockwell  International  Science  Center 
Thousand  Oaks,  California  91360 


ABSTRACT 


The  crack  closure  phenomenon  In  fatigue  under  plane  strain  is  studied 
by  means  of  an  inclined  strip  yield  superdislocation  model  and  finite  element 
stress  analysis,  with  particular  emphasis  on  the  effect  of  plastic  zone  con¬ 
strain  on  crack  closure.  Both  analyses  show  that  the  crack  closure  stress 
(normalized  with  respect  to  the  maximum  applied  stress)  increases  almost 
linearly  with  i/a  at  low  i/a  values  and  rises  more  slowly  with  increasing  z/a 
ratio  at  larger  z/a  values  ( l  being  slip  band  length,  a  the  half  crack  width). 


I.  INTRODUCTION 


Since  Elber's  [1,2]  discovery  of  the  crack  closure  and  the  introduction 
of  the  effective  stress- intensity  range,  a =*  Kmax  -  Kop,  where  Kfflax  Is  the 
maximum  stress-intensity  factor  and  Kop  Is  the  stress-intensity  factor  at  which 
the  crack  opens,  there  have  been  many  studies,  both  experimental  and  theoreti¬ 
cal,  on  crack  closure  and  its  effect  on  fatigue  crack  growth.  While  some 
Investigators  [3]  provided  experimental  evidence  that  crack  closure  is  a  surface 
effect,  others  [4-6]  concluded  that  closure  Is  a  general  phenomenon.  There  also 
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have  been  studies  [7,8]  which  Indicate  that  crack  closure  can  be  quite  pro¬ 
nounced  in  the  threshold  region  where  the  shear  fracture  mode  is  prevalent,  but 
the  arguments  presented  are  not  conclusive.  Comprehensive  theoretical  analyses 
of  closure  in  fatigue  crack  growth  based  on  the  Dugdale  model  in  plane  stress 
case  have  been  reported  recently  [9-11].  There  is  a  need  to  perform  similar 
theoretical  analyses  under  plane  strain  conditions.  This  is  one  of  the  prin¬ 
cipal  purposes  of  the  present  study.  Numerical  analyses  of  crack  closure 
employing  both  the  Inclined  strip  yield  superdislocation  model  [12-14]  and  a 
finite  element  method  were  used  in  this  investigation.  Another  objective  of 
this  study  concerns  the  crack  closure  behavior  in  polycrystalline  materials 
where  an  advancing  microcrack  runs  toward  a  grain  boundary  or  similar  obstacle 
which  restricts  the  slip  band  length  and/or  plastic  zone  size  to  the  dimension 
between  the  advancing  crack  tip  and  the  obstacle.  Literature  is  rich  in 
experimental  information  concerning  the  effect  of  grain  size  on  fatigue  crack 
closure  [15,16].  A  theoretical  understanding  of  the  experimental  observation  on 
grain  size  effects  is  important  to  the  kinetics  of  fatigue  crack  nucleation  and 
growth  in  real  materials. 


II.  INCLINED  STRIP  YIELD  SUPERDISLOCATION  ANALYSES 


The  Inhomogeneous 
advantageously  represented 
the  maximum  shear  stresses 
angle  to  the  crack.  811  by 


plastic  strain  produced  at  a  crack  tip  can  be 
on  a  microscale  by  dislocation  arrangements,  where 
are  symmetrically  located  on  planes  inclined  at  some 
and  Swinden  [17]  first  reported  the  results  of 
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numerical  calculations  for  dislocations  inclined  at  45°  to  the  crack  plane. 

Vitek  [18]  and  Riedel  [19]  have  recently  obtained  more  comprehensive  results  for 
slip  planes  Inclined  at  various  angles.  The  basic  equations  for  the  representa¬ 
tion  of  crack  tip  plasticity  based  on  an  inclined  strip  yield  superdislocation 
model  have  been  derived  by  Atkinson  and  Kanninen  [12,13]  and  used  by  Evans 
[14].  The  complex  variable  approach  of  Muskhelishvili  [20]  was  adopted  by  these 
investigators.  For  the  crack  and  slip  superdislocation  configuration  shown  in 
Fig.  1,  the  crack  of  width  2a  lies  along  the  X-axis  with  the  tips  of  the  crack 
situated  at  the  positions  (a,o)  and  (-a,o),  respectively.  The  superdislocation 
n,  where  n  *  1  to  m,  situated  at  positions  Zn  =  a  +  in  exp (la)  with  Burgers 
vectors  Bn,  experiences  a  shear  stress  -rn  given  by  each  of  the  m  discrete 
superdislocations. 


■  SHn  *  TO 


I, Bj <Gn j  *V'  n  ■  1,2  . . .m 


S  Is  the  applied  tensile  stress  normal  to  the  crack  plane,  v  is  the  shear 
modulus,  y  is  the  Poisson's  ratio,  Bj  is  the  Burger's  vector  of  the  jtfl  super- 
dislocation  and  the  remaining  terms  are  defined  below  (see  Eqs.  (7),  (8)  and 
(9)).  It  is  assumed  that  the  material  has  an  intrinsic  resistance  to  superdis¬ 
location  motion  along  the  slip  plane,  denoted  by  (the  friction  stress).  The 
equations  of  force  equilibrium  that  form  the  basis  of  the  mathematical  model  are 
formed  by  replacing  tn  of  Eq.  (1)  by  for  each  of  the  m  superdislocations. 

In  general,  the  interaction  between  the  crack,  applied  stress,  and  the 
superdislocations  produces  a  stress  field  behaving  like  K near  the  crack 


•.  j  *  • 
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tip,  where  r  is  the  distance  from  the  crack  tip  and  Kj  is  the  mode  I  stress- 
intensity  factor  given  by, 


Kj  *  S/^a  -  Kf  (2) 

is  the  component  of  the  stress-intensity  factor  that  results  from  the  effect 
of  the  stress  fields  of  the  superdislocations  acting  on  the  crack. 

Kj - -  1  B.F.  (3) 

4/»a(l  -  y)  j*l  J  J 

where  Fj  is  defined  below  (see  Eq.  (10))  and  the  other  terms  are  as  defined 
previously.  The  component  of  stress  containing  the  square-root  singularity  can 
be  suppressed  by  setting 

m 

S  '  4,.lf  -  rl  ^  ,4> 

In  our  representation  of  the  crack  tip  plasticity,  only  a  single 
superdislocation  rather  than  m  superdislocations  was  used.  This  representation 
is  justified  for  the  following  reasons.  First,  extensive  calculations  employing 
m  superdislocations,  where  m  is  a  variable  integer  used  in  the  manner  of 
Atkinson  and  Kannlnen  [12,13]  and  Evans  [14],  failed  to  yield  any  more  under¬ 
standing  than  that  according  to  a  single  superdislocation  representation. 
Furthermore,  numerical  solution  of  the  force  equations  becomes  considerably  more 
simplified  with  the  single  superdislocation  representation.  In  the  single 
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superdislocation  representation,  m  =  1,  Eqs.  (1)  and  (4)  are  reduced  to 


»  SH  + 


3'ttTT'-- 7T 


B  (G  +  K) 


(5) 


$  ■  T^rnr rT  BF 


(6) 


The  parameters  H,  G,  K  and  F  in  Eqs.  (5)  and  (6)  for  a  single  superdislocation 
with  Burger's  vector  B  situated  at  Z  =  a  +  i  exp( i 8)  are  given  by  [12], 


H  =  Sin  9 


I 

)  2 
(cos  9  +  a  z  Re 

e2l9 

2  T377 

JZ2  -  a2) 

(7) 


.21 0 


G  *  2i  sin  9  Re 


(Z2  -  a2) 


177 


a/2 


2Z2  +  a2 


Z  (l2  -  a2) _ 

(Z2  -  a 2 ) ( Z  -  Z)  2(z2  _  a2)3/Z 


.  (z2  -  An  -fz2  -  a2) 

.  «a219  | 

r  a2  Z  e19 

(Z  -  1)Z 

_ 

1  (z2  -  a2) 
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(9) 


Re 

ft*  -  a>)1/2l 

+  at  Re 

e19 

z  -  a 

71 — 7W? 

[ 

_(Z  -  a) (z2  -  ar) 

, 

(10) 


where  Re  is  the  real  part  of  and  1  Is  the  complex  conjugate  of  Z. 

For  small  scale  yielding,  i  «  a,  the  above  equations  are  further 
simplified.  Atkinson  and  Kanninen  [12]  and  Evans  [14]  treated  their  cases  in 
small  scale  yielding  approximation.  However,  for  slip  band  length  %  of  the 
order  of  half  crack  width  a,  small  scale  yielding  approximation  is  no  longer 
applicable  and  the  full  expressions  given  by  Eqs  (7)  to  (10)  must  be  used. 

When  the  applied  stress  S  Is  reduced  by  AS  to  (S  -  AS),  there  is  a 
tendency  for  the  dislocations  to  move  back  toward  the  crack,  but  the  movement 
will  be  resisted  by  the  friction  stress.  Only  those  dislocations  lying  between 
i  *  0  and  t  *  ir  (where  ir  Is  the  reversed  plastic  zone  dimension)  will  actually 
move.  When  a  new  equilibrium  Is  established,  the  dislocations  in  the  reversed 
plastic  zone  will  experience  a  friction  stress  equal  to  -t^.  Within  the  region 
0  <  i  <  ir,  the  net  strength  of  the  dislocations  becomes  (B  -  Br),  and  Eq.  (5) 
becomes 

-Tl  *  (S  -  A$)H  +  4ir-(Iu,  yy  (B  -  Br)  (G  +  K)  (11) 
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Taking  the  difference  between  Eq.  (5)  and  Eq.  (11),  one  obtains 

2tt  =  ASH  +  4V(T^  -y1  8r  (G  ♦  K)  (12) 

where  Br  can  be  viewed  as  the  contribution  to  dislocation  distribution  from  the 
reversed  shear  as  a  result  of  the  lowering  of  applied  stress  from  S  to  (S  - 
AS).  When  AS  3  S,  or  (S  -  AS)  *  0,  Br  from  Eq.  (12)  yields  directly  the  extent 
of  crack  closure  at  zero  applied  stress.  Combination  of  Eqs.  (5),  (6)  and  (12) 
yields  the  crack  closure  stress  Scc.  The  stress  Parameters  t. ,  S,  aS  and  the 
shear  modulus  v  In  Eqs.  (5)  to  (12)  are  made  dimensionless  by  expressing  these 
parameters  In  units  of  .  These  equations  could  not  be  solved  analytically. 

The  Interesting  parameters  such  as  B,  B  ,  S  ,  S  for  given  Inputs  of  crack  width 
(2a),  slip  band  length  U),  friction  stress  (t^),  slip  plane  angle  (a),  and 
shear  modulus  (u)  can  be  readily  obtained  by  numerical  iteration.  Typical  plots 
of  B,  8r,  versus  S/t^  for  aluminum  alloys  are  shown  In  Figs.  2  and  3, 
respectively,  for  9  =  70°  and  2a  *  100,  150  microns.  The  corresponding  solu¬ 
tions  of  z,  satisfying  the  singularity  suppression  Eq.  (6)  for  the  respective 
values  of  B  and  S/ij,  are  also  shown  In  Figs.  2  and  3  as  dotted  curves. 

In  order  to  evolve  the  change  of  crack  closure  stress  Scc/Smax 
(normalized  with  respect  to  maximum  applied  stress  Smax)  with  slip  band  lengr< 

1/ a  (normalized  with  respect  to  half  crack  width  a)  from  the  plots  shown  In 
Figs.  2  and  3,  It  Is  necessary  to  make  some  assunptions  In  regard  to  the  change 
of  friction  stress  with  slip  band  dislocation  density.  The  slip  band 
thickness  Is  assumed  to  Increase  In  direct  proportion  to  its  width,  so  that  for 
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a  given  superdislocation  Burger's  vector  the  average  dislocation  density  Is 
inversely  proportional  to  the  square  root  of  the  slip  band  length.  If  one 
further  assumes  that  the  friction  stress  to  superdislocation  motion  is  directly 
proportional  to  the  slip  band  dislocation  density  [21],  It  follows 


Tj  /  Tg  »  /lB/t)1/(B/i)0  (13) 

where  t0  and  (B/i)0  are  some  arbitrarily  chosen  reference  state  corresponding  to 
friction  stress  t0  and  slip  band  width  iQ.  The  crack  closure  stresses  calculated 
according  to  Eq.  (13)  are  plotted  versus  the  i/a  values  for  8  3  70°  as  solid 
curve  In  Fig.  4.  The  results  from  finite  element  analyses  (see  next  section) 
are  also  shown  in  Fig.  4  as  a  dotted  curve.  Experimental  data  [16]  for  a  high 
strength  aluminum  alloy  from  our  laboratory  are  presented  in  Fig.  4  as  solid 
circles. 


III.  FINITE  ELEMENT  ANALYSES 


The  NASTRAN  [22]  cracked  plate  piecewise  linear  analysis  program  with 
triangular  (CTRMEM)  elements  was  used  In  this  study.  The  cracked  plate  has  two 
planes  of  symmetry,  so  only  one  quarter  of  the  plate  needs  to  be  structured  and 
analyzed.  A  schematic  diagram  of  the  elements  representing  one  quarter  of  a 
cracked  plate  Is  shown  In  Fig.  5.  The  total  number  of  nodal  points  and  tri¬ 
angular  elements  used  varies  with  the  plate  width  (In  the  X-directlon)  to  the 
crack  width  ratio  studied.  The  right  edge  of  the  cracked  plate  is  constrained 
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from  moving  In  the  X-dlrectlon  to  simulate  a  fixed  boundary.  The  analysis 
Involves  stressing  the  plate  In  the  Y-dlrectlon  In  Increments  and  recalculating 
the  material  properties  for  each  element  as  a  function  of  the  element  stress  for 
the  last  load  Increment.  Plastic,  rather  than  nonlinear  elastic,  behavior  Is 
assuned.  The  theoretical  basis  of  two  dimensional  plastic  deformation  used  in 
NASTRAN  Is  that  developed  by  Swedlow  [23].  A  typical  uniaxial  stress-strain 
curve  for  7075-T6  aluminum  alloy  was  used  In  the  present  investigation.  Upon 
Incremental  loading  to  the  maximum  stress  S,,^,  loading  was  decreased  Incre¬ 
mentally  until  the  first  contact  of  the  cracked  nodal  point  or  points,  at  which 
the  applied  stress  is  called  the  crack  closure  stress  Scc.  Results  for  Scc/Smax 
are  plotted  versus  i/a  values  (where  i  *  w  -  a,  w  being  the  half  plate  width, 
and  a  the  half  crack  width)  In  Fig.  4  as  dotted  curves  for  comparison  with  the 
solid  curves  from  Section  II  according  to  the  superdislocation  model.  In  the 
finite  element  analyses,  Smax  was  chosen  to  be  nine-tenth  of  the  yield  stress  of 
the  material  In  all  the  calculations. 


IV.  DISCUSSION 


In  Figures  2  and  3  the  ratio  of  the  applied  stress  to  friction  stress 
Is  plotted,  regardless  of  the  state  of  the  friction  stress.  During  cyclic 
stressing,  the  friction  stress  during  the  loading  cycle  may  differ  from  that 
during  the  unloading  cycle.  The  effect  of  variation  of  the  friction  stress  tj, 
can  thus  be  readily  taken  Into  consideration  by  the  appropriate  modification  of 
the  values.  Although  these  calculations  are  based  on  a  very  simplified 
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superdislocation  model,  they  do  provide  theoretical  support  of  the  existence  of 
crack-closure  phenomenon  In  fatigue  crack  growth  under  plane  strain  condi¬ 
tions.  It  should  be  noted  In  Figs.  2  and  3  that,  for  each  data  point  of  the 
superdislocation  Burger's  vector  versus  S/t^  curve,  there  Is  a  unique  slip  band 
width  which  satisfies  the  singularity  suppression  Eq.  (6).  The  corresponding 
slip  band  width  values  are  plotted  versus  S/t^  also  In  Figs.  2  and  3  as  dotted 
curves.  The  predicted  dependence  of  crack  closure  on  slip  band  width,  based  on 
a  simplified  friction  stress  model  according  to  Eq.  (13)  (see  Fig.  4),  is  In 
reasonably  good  agreement  with  the  finite  element  analyses  results  described  in 
Section  III  and  with  experimental  observations  [16]. 
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Schematic  diagram  of  crack  and  superdislocation  configuration  (the 
inclined  strip  yield  zone  model). 


Superdislocation  Burger's  vector  (B)  and  slip  band  length  U/a, 
normalized  with  respect  to  half  crack  width)  versus  S/t<  (S  is  applied 
stress,  Tf  is  friction  stress)  for  a  *  50  microns,  8  *  /0#. 


Same  as  Fig.  2,  a  *  75  microns,  8s  70°. 


Crack  closure  stress,  Sc C/Sm-X,  normalized  with  respect  to  maximum 
applied  stress,  versus  slip  Band  length,  z/a,  normalized  with  respect 
to  half  crack  width  a,  for  8  *  70*. 


Nodal  element  configuration  for  finite  element  analyses. 
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